Laser field absorption in self-generated electron-positron pair plasma 
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Recently much attention has being attracted to the problem of limitations on the attainable 
intensity of high power lasers [A.M. Fedotov et al. Phys. Rev. Lett. 105, 080402 (2010)]. The laser 
energy can be absorbed by electron-positron pair plasma produced from a seed by strong laser field 
via development of the electromagnetic cascades. The numerical model for self-consistent study of 
electron-positron pair plasma dynamics is developed. Strong absorption of the laser energy in self- 
generated overdense electron-positron pair plasma is demonstrated. It is shown that the absorption 
becomes important for not extremely high laser intensity / ~ 10 24 W/cm 2 achievable in the nearest 
future. 
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Due to an impressive progress in laser technology, laser 
pulses with peak intensity of nearly 2 x 10 22 W/cm 2 are 
now available in the laboratory When the matter 
is irradiated by so intense laser pulses ultrarelativistic 
dense plasma can be produced. Besides of fundamen- 
tal interest, such plasma is an efficient source of particles 
and radiation with extreme parameters that opens bright 
perspectives in development of advanced particle acceler- 
ators , next generation of radiation sources 0, 0] , labo- 
ratory modeling of astrophysics phenomena 0], etc. Even 
higher laser intensities can be achieved with the coming 
large laser facilities like ELI (Extreme Light Infrastruc- 
ture) @ or HiPER (High Power laser Energy Research 
facility) Q . At such intensity the radiation reaction and 
quantum electrodynamics (QED) effects become impor- 
tant [l-[l3|. 

One of the QED effects, which has recently attracted 
much attention, is the electron-positron pair plasma 
(EPPP) creation in a strong laser field [III El- The 
plasma can be produced via avalanche- like electromag- 
netic cascades: the seed charged particles are acceler- 
ated in the laser field, then they emit energetic pho- 
tons, the photons by turn decay in the laser field and 
create electron-positron pairs. The arising electrons and 
positrons are accelerated in the laser field and produce 
new generation of the photons and pairs. It is predicted 
12] that an essential part of the laser energy is spent on 
EPPP production and heating. This can limit the attain- 
able intensity of high power lasers. That prediction was 
derived using simple estimates, therefore self-consistent 
treatment based on the first principles is needed. 

The collective dynamics of EPPP in strong laser field 
is a very complex phenomenon and numerical modeling 
becomes important to explore EPPP. Up to now the nu- 
merical models for collective QED effects in strong laser 



field have been not self-consistent. One approach in nu- 
merical modeling is focused on plasma dynamics and ne- 
glects the QED processes like pair production in the laser 
field. It is typically based on particle-in-cell (PIC) meth- 
ods and uses equation for particle motion with radia- 
tion reaction forces taken into account The second 
one is based on Monte Carlo (MC) algorithm for photon 
emission and electron-positron pair production. This ap- 
proach has been used to study the dynamics of electro- 
magnetic cascades jl4|. However, it completely ignores 
the self-generated fields of EPPP and the reverse effect of 
EPPP on the external field. The latter effect is especially 
important to determine the limitations on the intensity 
of high power lasers 12 , 1|| . 



Quantum effects in strong electromagnetic fields can 
be characterized by the dimensionless invariants la, [17 1 
Xe = ehj{m i c^)\F llv 'p v \ w j(F±/eE cr ) and Xi ~ 
(fku/mc 2 )(F±/eE cr ), where is the field-strength ten- 
sor, is the particle four-momentum, fnu is the pho- 
ton energy, 7 is the electron gamma-factor, F± is the 
component of Lorentz force, which is perpendicular to 
the electron velocity, E cr — m 2 c 3 /(eh) = 10 16 V/cm is 
the so-called QED characteristic field, ft is the Planck 
constant. Xe determines photon emission by relativistic 
electron while x 7 determines interaction of hard photons 
with electromagnetic field. QED effects are important 
when Xe ^ 1 or x 7 > 1. If % e > 1 then hui ~ 7171c 2 and 
the quantum recoil imposed on the electron by the emit- 
ted photon is strong. The probability rate of emission 
of a photon with energy Huj by relativistic electron with 
gamma-factor 7 can be written in the form 17- 19j] 
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where hu> is the photon energy, m is the electron mass, 
c is the speed of light, S = 2£/[3(l — £)Xe] and £ = 
hjj/(^jmc 2 ). hjjdW em can be considered as the energy 
distribution of the electron radiation power. For electron 
radiation in constant magnetic field B perpendicular to 
the electron velocity it reduces to the synchrotron radia- 
tion spectrum in the classical limit Xe "C 1 0, [2(| ■ The 
probability rate of electron-positron pair pr oduction by 
decay of a photon with energy hui is [17Hl9j 
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where 6 = 2/ (3x-yV- r l+)i V- = imc 2 / (fuv) and r/ + = 

1 — 77 are the normalized electron and positron energies, 

respectively. It follows from Eq. (|2]) that in the classical 
limit Xi ^ 1 this probability is exponentially small. 

To study EPPP dynamics we have developed two- 
dimensional numerical model based on PIC/MC meth- 
ods. Similar methods has been used previously for 
modeling of discharges in gases 21]. Recently a one- 
dimensional PIC/MC model has been developed to sim- 
ulate pair cascades in magnetosphere of neutron stars 
22]. However the latter model is electrostatic and the 
classical approach is used for photon emission with ra- 
diation reaction force in the equation of motion. In our 
numerical model we use more general approach. We ex- 
ploit the fact that there is a large difference between the 
photon energy scales in EPPP. The photon energy of the 
laser and plasma fields is low (Hu> <C mc 2 ) while the en- 
ergy of the photons emitted by accelerated electrons and 
positrons is very high (hu> ^> mc 2 ). The emitted photons 
are hard and can be treated as particles. Conversely, the 
evolution of the laser and plasma fields is calculated by 
numerical solution of the Maxwell equations. Therefore, 
the dynamics of electrons, positrons and hard photons as 
well as the evolution of the plasma and laser fields are cal- 
culated by PIC technique while emission of hard photons 
and pair production are calculated by MC method. 

The photon emission is modeled as follows. On every 
time step for each electron and positron we sample a pho- 
ton emission by a probability distribution which approx- 
imates Eq. (fTJ) with the accuracy within 5%. The new 
photon is included in the simulation region. The coordi- 
nates of a new photon are equal to the electron (positron) 
coordinates at the emission instance. The photon mo- 
mentum is parallel to the electron (positron) momentum. 
The electron (positron) momentum value is decreased by 
the value of the photon momentum. Similar algorithm 
is used for modeling of pair production by photons. The 
new electron and positron are added in the simulation 
region while the photon that produced a pair is removed. 



The sum of the electron and positron energy is equal to 
the photon energy. The pair velocity is directed along 
the photon velocity at the instance of creation. 

The MC part of our numerical model has been bench- 
marked to simulations performed by other MC codes. 
We simulated the electromagnetic showers in a static ho- 
mogeneous magnetic field, the interaction of relativistic 
electron beam with a strong laser pulse, and the devel- 
opment of electromagnetic cascades in circularly polar- 
ized laser pulses. The obtained results are in reasonably 
good agreement with those published by other authors 
and are discussed in Ref. |19| . The particle motion and 
evolution of the electromagnetic field are calculated with 
standard PIC technique ,23;] . The PIC part of the model 
is two-dimensional version of the model used in Ref. [2~3 | . 
In order to prevent memory overflow during simulation 
because of the exponential growth of particle number 
in a cascade, the method of particles merging is used 
[22| . If the number of the particles becomes too large the 
randomly selected particles are deleted while the charge, 
mass, and the energy of the rest particles are increased 
by the charge, mass, and energy of the deleted particles, 
respectively. 

We use our numerical model to study production 
and dynamics of EPPP in the field of two colliding 
linearly polarized laser pulses. The laser pulses have 
Gaussian envelopes and propagate along the x-axis. The 
components of the laser field at t = are 
a exp(-y 2 /cr 2 )sinC e -0+*o) 2 /°-2 ± e -{x-x f/al 

where the field strengths are normalized to mcu>L/ |e|, 
the coordinates are normalized to c/u>l, time is nor- 
malized to 1/ojl, o-q — \e\Eo/(mcu>L,), Eq is the electric 
field amplitude of a single laser pulse, is the laser 
pulse cyclic frequency, 2xq is the initial distance be- 
tween the laser pulses, £ = x — tf>, and <f> is the phase 
shift. The parameters of our simulations arc <p = 0.8-7T, 
do = 1.2 • 10 3 , a x = 125, oy = 40, x = a x /2 that for 
the wavelength A = 2itc/ujl = 0.8 /im correspond to the 
intensity 3 • 10 24 W/cm 2 , pulse duration 100 fs, the focal 
spot size 10 /Ltm at 1/e 2 intensity level. The cascade is 
initiated by a single electron located at x = y = with 
zero initial momentum for t = when the laser pulses 
approach to each other (the distance between the pulse 
centers is a x ). 

The later stage (t = 25.5A/c) of the cascade develop- 
ment is shown in Fig. [1] where the electron and photon 
density distributions and the laser intensity distribution 
are presented. The laser pulses passed through each other 
by this time instance and the distance between the pulse 
centers becomes about 1.6a x . It is seen from Figs. Q] that 
the micron-size cluster of overdense EPPP is produced 
and the laser energy at the backs of the incident laser 
pulses is spent on EPPP production and heating. The 
plasma density exceeds the relativistic critical density 
aQn cr in about 2 times, where n cr — muj I (87re 2 ) is the 



E y , B z 



3 



1 r~ 



-i 1 r 



o - 



-5 - 



- 



-5 - 



_i I i_ 

-i 1 13 



_i I i_ 

r — 1 — r 



_i I 1 I 1 L 



pi 



5» 



5 - 

- 



-5 - 



~i 1 r 



~> 1 r 



_l I I I I I L. 



Pi 



-40 -20 



20 40 



x/X 



FIG. 1: The normalized electron density p e = n e /(aon cr ) 
(a), the normalized photon density p 7 = n 7 /(aon cr ) (b) and 
the laser intensity normalized to the maximum of the initial 
intensity pi (c) during the collision of two linearly-polarized 
laser pulses at t — 25.5A/c. 




FIG. 2: The electron and positron energy (solid line), the pho- 
ton energy (dotted line), the laser energy (dashed line) and 
the total energy of the system (dash-dotted line) as functions 
of time. All the energies are normalized to the initial energy 
of the system. 



nonrelativistic critical density for the electron-positron 
plasma. The evolution of the particle and laser energy 
is shown in Fig. [2] It is seen from Fig. [2] that about 
a half of the laser energy is absorbed by self-generated 
EPPP and then mostly reradiated in ultrashort pulse of 
gamma-quanta. The total energy of the particles in the 




FIG. 3: The number of the electrons produced in the cascade 
(line 1) and the EPPP density normalized to the relativistic 
critical density (line 2) as functions of time. 




FIG. 4: The profile of the electron density along x-axis at 
y = for initial stage ct = 6.4A (line 1) and for the later 
stage ct — 16. 6A (line 2) of the cascade development. The 
electron density for ct = 16. 6A is normalized to aon cr and 
that for ct = 6.4A is normalized to 3.3 x W~ 6 aon cr . 



cascade and the electromagnetic field is conserved with 
accuracy about 1% during our simulation. 

At initial stage of the cascade development the num 
ber of created particles is growing exponentially N ~ e 
(l2| . where T is the multiplication rate. It follows from 
the energy conservation law that the number of parti- 
cles that can be created is limited by the laser pulses 
energy. Thus, at some instant the exponential growth 
is replaced by much slower growth. Equating the initial 
energy of laser pulses to the overall particles energy af- 
ter the pulse collision we get N ~ aQa x a^N /j, where 
we assume N e ~ N p ~ N p h ~ N, N e , N p and N p h are 
the number of electrons, positrons and photons produced 
by the cascade, respectively, mc 2 7 is the average parti- 
cle energy, = n cr (c/oj) 3 = A/(167rr e ), r e — e 2 /(mc 2 ). 
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FIG. 5: The dependence of Xe for the primary electron on 
time. 



The multiplication rate decreases when the field strength 
goes down, that, by turn, occurs if the plasma density 
reaches the value aon cr . This is in good agreement with 
the numerical results shown in Fig. [3J where the mul- 
tiplication rate drops dramatically and EPPP density 
reaches the value about aon cr at the same instant of 
time t s f=a lOA/c. The value of t s can be estimated as 
t s w r-MniV. It follows from Fig. U that T w 0.6w L 
for t < t s . The typical lifetime t em for electrons and 
positrons with respect to hard photon emission can be 
estimated as 1/r > 1/ujl [HI- Thus, for the parameters 
of numerical simulation 7 can be estimated as 7 ~ ao 
hence N ~ aoa x a^N ~ 4 • 10 14 and ct s /X ~ 9 that are 
in good agreement with the corresponding values from 
Fig.H 

It is shown in Ref. [HI that the cascades do not arise 
in B-nodc of linearly polarized standing electromagnetic 
wave so far as Xe and Xi are l ess than unity. However, 
our numerical simulations show that the cascade quasi- 
periodically develops between B and E nodes of such a 
wave. This is because under such conditions the electron 
motion becomes complicated and is not confined to the 
direction of polarization on the temporal scales about 
the laser period. It turns out that there occur the time 
intervals of duration of the order of u;^ 1 with % e > 1 [17 
(see Fig. [5]) on which the cascade can develop. The time 
modulations of N e (t) and of EPPP density along the x- 
axis at the initial stage of the cascade development are 
seen in Fig. [3] and 21 respectively. At the later stages the 
spatial modulation of the density is strongly smoothed 
out due to EPPP expansion (see Fig. [4] line 2). 

In conclusion we develop the numerical model which al- 
lows us to study EPPP dynamics in strong laser field self- 
consistently. We have demonstrated efficient production 
of EPPP at the cost of the energy of the laser pulses. We 
show that even not extremely high intensity laser pulses 
(/ ~ 10 24 W/cm 2 with duration ~ 100 fs) can produce 



overdense EPPP so that the QED effects can be experi- 
mentally studied with near coming laser facilities like ELI 
Q and HiPER Q. The simulations and estimates show 
that for intensity I > 10 26 W/cm 2 the overdense EPPP 
can be produced during a single laser period. In such 
high-intensity regime few-cycle laser pulses can be used in 
experiments. High-energy photons or electron-positron 
pair can be also used as a seed to initiate cascade instead 
of an electron. Photon-initiated cascade can be more 
suitable for experimental study in low intensity regime 
(/ ~ 10 24 W/cm 2 ) because the laser intensity threshold 
for pair creation in vacuum is about ~ 10 26 W/cm 2 [25| . 
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